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An assessment of an acoustic analogy for the mixing noise component of jet noise in the 
presence of an infinite surface is presented. The reflection of jet noise by the ground changes 
the distribution of acoustic energy and is characterized by constructive and destructive 
interference patterns. The equivalent sources are modeled based on the two-point cross- 
correlation of the turbulent velocity fluctuations and a steady Reynolds- Averaged Navier- 
Stokes (RANS) solution. Propagation effects, due to reflection by the surface and refraction 
by the jet shear layer, are taken into account by calculating the vector Green’s function of 
the linearized Euler equations (LEE). The vector Green’s function of the LEE is written in 
relation to Lilley’s equation; that is, approximated with matched asymptotic solutions and 
the Green’s function of the convective Helmholtz equation. The Green’s function of the 
convective Helmholtz equation for an infinite flat plane with impedance is the Weyl-van 
der Pol equation. Predictions are compared with an unheated Mach 0.95 jet produced by 
a nozzle with an exit diameter of 0.3302 meters. Microphones are placed at various heights 
and distances from the nozzle exit in the peak jet noise direction above an acoustically 
hard and an asphalt surface. The predictions are shown to accurately capture jet noise 
ground effects that are characterized by constructive and destructive interference patterns 
in the mid- and far-field and capture overall trends in the near-field. 
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Nomenclature 

Constant associated with dilatation source 

Constant associated with unsteady force per unit volume source 
Speed of sound 

Constant of integral length scale 
Constant of integral velocity scale 
Constant of integral time scale 
Nozzle exit diameter 

Variation of two-point cross-correlation of velocity fluctuations with space and time 

Rudwick boundary loss factor 

Frequency 

Unsteady force per unit volume associated with velocity fluctuations 

Green’s function of the convective Helmholtz equation 

Green’s function of Lilley’s equation 

Location of turning point 

Turbulent kinetic energy 

Integral length scale in the streamwise direction 
Integral length scale in the y direction 
Integral length scale in the z direction 
Mach number 

Fully expanded Mach number 
Pressure 

Reflection coefficient 
Flow resistivity 
Radial direction 
Spectral density 
Total temperature ratio 
Time 
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Streamwise velocity component in the x direction 

Integral scale of velocity 

Radial velocity component in the y direction 

Radial velocity component in the 2 : direction or numerical distance 
Observer position 
Streamwise direction 
Source position 

Vector from source to observer 
Surface impedance 
Ratio of specific heats 
Dirac delta function 

Dissipation rate of turbulent kinetic energy 
Vector between two source locations (£ , y, C) 

Dilatation rate or observer angle from nozzle downstream axis 
Logarithm of the pressure 

Vector Green’s function of the linearized Euler equations 

Density 

Retarded time 

Integral time scale 

Azimuthal angle 

Observer angle from the nozzle inlet axis 
Specific dissipation rate of turbulent kinetic energy 
Radian frequency 


Jet fully expanded quantity 
Property of the primary jet 
Property of turbulence 
Ambient value 


Introduction 

An acoustic analogy is developed based on the Euler equations for the prediction of jet mixing noise in the 
vicinity of an arbitrary geometry. The analogy is based on the development by Morris and Farassat 1 but has 
unique differences. The Euler equations, written in terms of the logarithm of the pressure, are rearranged into 
a left hand side operator consistent with the linearized Euler equations (LEE) and right hand side equivalent 
sources. The far-held pressure is written in terms of an integral solution of the governing equations. It 
appears as a volumetric integral of the equivalent sources and the vector Green’s function of the LEE. 
Any adjoint vector Green’s function solver for the LEE can be used within the analogy. Equivalent source 
models are created for both the dilatation and unsteady force per unit volume source terms. The equivalent 
sources follow the models of Tam and Auriault 2 & Morris and Boluriaan . 3 The model is dependent on 
steady Reynolds-Averaged Navier-Stokes (RANS) solutions but could easily be based on an empirical how- 
helcl or unsteady simulation. This allows the nozzle pressure ratio, temperature ratio, and nozzle geometry 
to be connected directly to the aerodynamic solution and resultant jet noise. A model of the two-point 
cross-correlation of the unsteady velocity fluctuations within the jet connects the turbulence statistics to the 
prediction. It is based on measurement and the turbulence statistics generated by a two-equation closure of 
the steady RANS equations. 

The refraction of the sound through the meanhow, diffraction by surfaces, or reflection by the ground 
or airframe, are all handled by the solution of the vector Green’s function of the LEE. For simplicity in 
this investigation, the vector Green’s function of the LEE is written in terms of Lilley’s 4 equation. An 
approximate solution of Lilley’s equation is formed about the Green’s function of the convective Helmholtz 
equation. Analytic forms of the Green’s function of the convective Helmholtz equation are used to account 
for sound reflection on an infinite flat surface which is representative of the ground. The impedance model 
of Delany and Bazley 5 is used to simulate the damping of the reflected acoustic waves, but can be neglected 
for an acoustically hard wall. The form of the Green’s function of the convective Helmholtz equation is the 
Weyl-van der Pol 6 formula. The geometry chosen to exercise the model coincides with the experiment by 
Miles . 7 

The methodology presented in this paper is developed to address many of the issues encountered by pre- 
vious investigators. Experimental measurements conducted near reflecting surfaces pose unique difficulties. 
These are summarized well by the studies of Seiner et al ., 8 Butzel , 9 Schlinker et al ., 10 Pao et al ., 11 Huber 
and Sogeti , 12 and many others, who performed experiments with microphones at or near various surface 
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locations relative to static jet and aircraft flyover. With careful microphone placement and corrections, 
accurate measurements can be conducted. 

A large body of prediction methodologies were developed to address reflection and impedance effects. In 
particular, Illston et a /. 13 performed a combined theoretical and experimental study on jet noise, erosion, 
and ground effects. Their methodology used a Harper-Bourne 14 like prediction approach by integrating 
the source strength per unit length along the centerline jet axis. The Aircraft Noise Prediction Program 
(ANOPP) (see Zorumski 15 for details), accounts for ground effects from various sources such as the fan, 
jet, and airframe. ANOPP uses ground reflection corrections described by Chien and Soroka 16 and the 
impedance model by Delany and Bazley . 5 Finally, McLaughlin et al. 1 ' performed a series of experiments 
and predictions using supersonic jets with and without a ground plane. They created a model to describe the 
source distribution within the jet plume and predictions showed that the interference patterns were correctly 
predicted at high frequencies. 

In the following sections, the model development is presented starting with the governing equations. An 
overview of the implementation of the method in a computer program, its connection to a steady RANS 
solver, a description of the model calibration, and method to find the steady RANS solution are shown. 
Results of the predictions are compared with the experiment of Miles.' 

Mathematical Model 

The acoustic analogy presented herein is based on that of Morris and Farassat . 1 The governing equations 
are the Euler equations, 
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where D / Dt is the material derivative, c is the speed of sound, u is the velocity, and 7 r = r y~ 1 ln{p/p 00 ). By 
linearizing the governing equations about a non-radiating base-flow and retaining terms on the left hand side 
that are linear in fluctuation then the inhomogeneous LEE are formed. 


and, 
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is a unsteady force per unit mass. In both the equivalent sources only the second order fluctuations are 
considered to contribute to the sound field. 

Green’s function of the LEE is defined as, 
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where D a /Dt = d/dt + Ujd/dxj. Let the vector Green’s function be periodic, 

/ OO 

■Kg{x\y; u>) exp [iu(t - r)]duj (9) 

-OO 

and 

ng(x\y,uj) = T J 7Tg(x, t\y,r) exp[—iuj(t — r)]dr. (10) 

The fluctuating far-field pressure is a convolution integral of the vector Green’s function and the equivalent 
sources, 

/ oo noo 3 

... / Tr° g (x,t\y,T)6(y,T) + Kg{x,t\y,T)f n (y,T)dTdy. (11) 

-OO J - OO n=1 

The spectral density is defined as the Fourier transform of the autocorrelation of the fluctuating pressure, 


S(x,u) = / p'(x, t)p'(x, t + t) exp[— iojT]d,T. 


Using Eqns. 11 and 12 yields, 
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By definition the periodic vector Green’s function of the LEE, 

Kg*(x,y,u) = n™(x,y,-u) 

and following Tam and Auriault, two closely placed source points in the jet are related by, 
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Using Eqns. 14 and 15 with Eqn. 13 yields, 
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We have now obtained a general formulation that involves the two-point cross-correlation of the equivalent 
sources and the vector Green’s function of the LEE. A model is required for the two-point cross-correlation 
of the equivalent sources and the vector Green’s function of the LEE. 

The dilatation rate source is modeled following Tam and Auriault 2 and written in the form of Morris 
and Boluriaan, 3 

0(y, T)0(y + y, t + t) = A 2 E(r), r) (17) 

' s 

where A s is a constant. The two-point cross-correlation of the unsteady force per unit volume is modeled 
following Morris and Boluriaan, 3 

fn(y,T)f m (y + y,t + T) = B ^ Us ^ c ^ Us (18) 
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where B s is a constant. A model must be formed for the space-time decay function E(t),t). Following 
Ribner’s 18 postulate, the spatial and temporal terms of E(rj, r) are separable, 


E(r], t) = exp 
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Equations 17 through 19 are used in Eqn. 16. Integrals involving the variables of integration r, y, £, and £ 
are performed analytically. Assumptions required for analytical integration involving r) and r require r s > 0, 
li > 0, Coo > 0, x > 0, and for these quantities to be real. These are reasonable assumptions physically as 
the turbulent kinetic energy at the boundary is zero and we restrict our modeling of the jet flow where x = 0 
represents the nozzle exit. These integrations result in, 
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The integral scales of turbulence in Eqn. 20, l x , l y , l z , u s , and r a , are required to make a prediction. 
These could be found by simple empirical models or unsteady CFD simulations. Here, they are related to a 
steady RANS solution by, 


lx(y) = ciK(y)/e(y), 


(21) 


r s (y) = c T K(y)/e(y), (22) 

and, 

u s (y) =c u y/2K(y)/3. (23) 

The coefficients are c T = 0.30, c u = 1.00, and q = 1.00. These coefficients are based on a reference jet 
operating at the sonic condition, TTR = 3.20, and a 0.0508 m convergent nozzle at the sideline location of 
R/ D = 100. They have been calibrated by the same methodology of Tam and Auriault. 2 Their variation from 
Tam and Auriault is attributed to the use of a different acoustic analogy and steady RANS solver. The cross- 
stream integral length scales, l y and l z , are set to 3/10 of l x , which corresponds to experimental observation. 
These coefficients are never altered irrespective of the nozzle Mach number, diameter, temperature, nozzle 
geometry, or any other parameter. 

A solution for 7r™ is now required. One method suitable for calculating the vector Green’s function of 
the LEE is by finding the Green’s function of Lilley’s 4 equation. The vector Green’s function of the LEE 
are related to the Green’s function of Lilley’s 4 equation by, 
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where the subscript of the source vector y denotes the direction the partial derivative is taken. The problem 
is now reduced to finding the Green’s function of Lilley’s equation instead of the full vector Green’s function 
of the LEE. Tam and Auriault 19 and Raizada 20 performed this technique with success. These methods used 
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an adjoint approach and are relatively computationally inexpensive. However, these methods do not include 
a scattering surface such as the airframe or ground, thus analytic or numerical approximations of Lilley’s 
equation that have the potential to include a scattering surface are sought. 

Analytical solutions of Lilley’s 4 equation can be found but only for very special mean flows. These 
are instructive but not helpful for more realistic jets. An asymptotic solution of Lilley’s equation for low 
frequencies was found by Goldstein 2 1,221 23 and for high frequencies by Balsa et al. 24 By examination of the 
asymptotic solutions, the Green’s function of the convective Helmholtz equation, g , and a matching function, 
an approximate Green’s function of Lilley’s equation is created, 
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9 

exp [-Cgw/Coo] g 


if Re[g* 0 \ > 0 
if Re[g*\ < 0 


(28) 


I (1 — u/Cao COS 9) 2 


(c 2 /cy 


— cos 2 9. 


9o = 

Co is Ribner’s convection coefficient, 

C 0 = \J{1 -u/c 00 cos9) 2 + cl a {u/c rx ) 2 

where Cq q = (4/25) a/ic,- /Cqo — 1. The convective amplification power coefficient is, 


(29) 


(30) 


c a = 3 — tanh [fDj/uj ] . (31) 

The benefit of this approach is existing numerical solvers can provide the Green’s function of the con- 
vective Helmholtz equation. For example, the Fast Scattering Code (FSC) by Tinetti and Dunn, 25 any 
Boundary Element Method (BEM) (see an overview by Katsikadelis 26 ), or FastBEM by Liu. 2 ' This permits 
the opportunity of finding the jet noise scattering from any surface as long as the aerodynamic noise sources 
are minimally affected by the scattering geometry. 

For static jet predictions the Green’s function of the Helmholtz equation is used, g = exp[*fc;z]/(47r;z). 
If an infinite plane exists parallel to the jet centerline then the tailored Green’s function of the Helmholtz 
equation is represented by the Weyl-van der Pol 6 formula, 
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where I? 2 is the distance from an equivalent source point on the opposite side of the infinite plane to the 
observer and R p is the reflection coefficient, 
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and F(w) is the Rudwick boundary loss factor, 


(33) 


F(w) = 1 + iy/nwex p [— w 2 ] erfc [— iw] (34) 

where erfc is the complementary error function. We numerically evaluate erfc using the identity erfc[— iw] = 
1 — erf[— iw}. w is the numerical distance, 

w=^(l + i)^/2TrfR 2 /c 00 [cos 9 + (pooCoo) -1 ] • (35) 

There are many models for the surface impedance Z s . Following ANOPP, 15 one appropriate choice is 
the model of Delany and Bazley, 5 

^ = PooCoo [1 + 0.0571A(/)-°- 754 + 70.87A(/)-°- 732 ] (36) 

where X(f) = Poof / R 3 - Rs is the flow resistivity. The asymptotic solution of Z s for R s approaching infinity 
yields the hard- wall solution. 
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Results 


The developed mathematical models are now assessed by comparing them with the experimental mea- 
surements of Miles. 7 Here, a convergent nozzle with a 0.3302 m (13 inch) exit diameter was placed in an open 
environment with the jet centerline axis 3.89 m (12.75 ft) above the ground. The nozzle centerline axis was 
parallel to the ground. The ground was described as asphalt and has a flow resistivity of R s = 27000 Pas/ 
m 2 , which is the average flow resistively described by Attenborough et al. 2S The ambient environment had a 
temperature of 288.9 K and was described as dry. Corrections have been made for atmospheric absorption. 
Microphones were placed relative to the nozzle exit at an angle ip = 145 degrees from the upstream axis. 
Three sets of microphones were located at this observer angle at ground shadow distances (distance from a 
point on the ground under the nozzle measured to a point under the microphone) of 7.62, 15.24, and 22.86 
m. At each of these observer distances, microphones were placed at 0.0061, 1.52, and 3.89 m above the 
ground. Lossless sound pressure level spectra were calculated in 1/3 octave bands from 50 to 20,000 Hz. 

The coordinate system shown in Fig. 1 illustrates the relative positions of the microphones with respect 
to the nozzle. The nozzle is shown on the left side of the figure and the jet plume develops from the left 
to right direction along the positive axaxis. The y - axis is normal to the ground and the ^-axis is parallel 
to the ground in the sideline observer direction. The height of the nozzle is specified by the distance along 
a normal from the ground plane to the center of the nozzle exit. The observer angle, ip, is the angle from 
the upstream nozzle centerline axis to the observer. The azimuthal angle, <p, is measured about the x-axis 
originating from the x-y plane. Values of x are specified to find the microphone position relative to the 
nozzle exit along with the microphone height. 

The Steady RANS Solution 

The developed prediction model requires the knowledge of a mean flow. The steady RANS solution is found 
with the Fully Unstructured Navier-Stokes (FUN3D) solver from The NASA Langley Research Center. For 
more information on FUN3D see Anderson and Bonhaus 29 or Nielsen. 30 The steady RANS equations are 
closed with the Menter 31 Shear Stress Transport (SST) turbulence model. Because the nozzle is axisymmetric 
and the ground is not specified as a boundary condition in the CFD calculation, the resultant steady flow- 
field is axisymmetric. The flow-field within the first quadrant of the y-z coordinates is obtained with 
the CFD calculation. After the calculation is complete, only the x-y plane is retained. The computer 
program that evaluates the prediction model rotates the steady RANS solution about the centerline axis 
when performing azimuthal integration. The methodology is compatible with other steady RANS solvers 
and has been demonstrated with the NPARC Alliance Wind-US solver (See Nelson 32 for details). 

The geometry of the nozzle is captured by the computational domain. The computational domain is 
unstructured and contains a region inside the nozzle, a region outside the nozzle upstream of the nozzle exit, 
and a region for the plume of the jet. The plume region extends 100 nozzle diameters downstream from the 
nozzle exit and 50 nozzle diameters in the radial direction from the jet centerline. A complete description 
of the methodology to construct the computational domains used for these types of simulations is found in 
Miller and Veltin. 33 The nozzle geometry included in the steady RANS simulation is based on the SMC000 
nozzle of NASA Glenn Research Center. See Bridges et al. 3i for information on the development of the SMC 
series of nozzles. The nozzle contour is shown in Fig. 2 along with a portion of the computational domain 
in the jet near-field. The number of grid points in the domain is 205,251 and the number of elements is 
204,000. 

Boundary conditions are chosen to replicate the experiment of Miles. 7 The 0.3302 m diameter convergent 
nozzle operates with a fully expanded Mach number of 0.95 and total temperature ratio of 1.00. This results 
in a fully expanded exit velocity of 301.66 m/s, fully expanded static temperature of 247.7 K, fully expanded 
speed of sound of 315.5 m/s, and fully expanded density of 1.43 kg/m 3 . The characteristic frequency of the 
jet is 913.6 Hz and the Reynolds number based on the fully expanded jet diameter is 8.944 x 10 6 . A portion 
of the jet flow-field predicted by FUN3D is shown in Fig. 3. Contours of the streamwise velocity component 
are shown in the top half plane and contours of turbulent kinetic energy are shown in the lower half plane. 
The spatially independent variables are normalized by the nozzle exit diameter. Doubling the number of 
grid points in the streamwise and cross-stream directions yield the same solution. Global residuals of the 
final CFD solution show reductions of over eight orders of magnitude. 
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Aeroacoustic Predictions 

A free-field calculation is required to validate the prediction method. Since experimental data are not 
available for this jet operating condition and nozzle diameter in the free-held, the result of the developed 
acoustic analogy is compared with the SAE ARP876 35 methodology. The SAE ARP876 35 prediction method 
for single stream shock-associated noise from convergent nozzles at supercritical conditions is included as a 
prediction module in NASA’s ANOPP. 15 The results are presented as lossless Sound Pressure Level (SPL) 
versus one-third octave band center frequency. In the sideline direction, the agreement between the prediction 
and SAE ARP876 is satisfactory at all frequencies as shown in Fig. 4. 

Predictions of jet noise about the ground plane are now presented. For these comparisons, the nozzle 
height is 3.89 m and the observer angle is if) = 145 deg. The ground shadow distance is defined as \Jx\ + x\ 
and is representative of the distance between the nozzle exit and microphone along the ground. The ground 
shadow distance between the nozzle exit and microphone is R/D = 69.23. Figure 5(a) shows the prediction 
of the jet noise compared with the experiment of Miles.' In this comparison the prediction of the jet noise 
above an acoustically hard surface (ideal) is denoted by a black line with circles, above an asphalt surface 
(impedance) as a blue line with deltas, and the measurement is denoted by a red line with squares. The 
y - axis is sound pressure level (SPL) in one-third octave bands and the x-axis is frequency. All following 
figures use this convention. The ground reflections create select destructive and constructive interference 
at various frequencies. For example, at 80 Hz the magnitude of the predicted and measured constructive 
interference agree well. An example of destructive interference occurs at 125 Hz. At frequencies higher than 
2500 Hz, the slopes of both the prediction and experiment are approximately the same. However, there is 
an under-prediction by about two dB through this range. 

A second comparison is shown in Fig. 5(b) at R/D = 69.23 but the microphone height is lowered to 
1.59 m. As the microphone is lowered towards the ground the constructive and destructive interference 
pattern increases in frequency. This trend matches the analytic solution of the Helmholtz equation with 
a point source above an infinite flat surface. Both predictions have matching constructive and destructive 
interference patterns, however these are not aligned with experiment. For example, the first measured 
destructive interference frequency is near 300 Hz but both predictions are located near 175 to 200 Hz. In 
addition, the amount of constructive interference in the prediction under-predicts that of experiment. The 
overall error between the hard-wall prediction and measurement is approximately two to four dB over 1000 
Hz and has considerable error approaching 10 dB at lower frequencies. This is consistent with the impedance 
ground surface prediction that is generally two dB lower over the same frequency ranges. 

When the observer is located near the ground, predictions become more challenging. The observer 
height is set at a height of 0.0061 m. The shadow distance between the microphone and nozzle remains at 
R/D = 69.23. A comparison between the predictions and experiment is shown in Fig. 5(c). The predictions 
show much more interference relative to experiment, which is probably due to the fact that microphones 
do not measure pressure at a point but over a small volume of space, and that no air gap exists between 
the microphone and the ground. The interference pattern is strongly affected by the microphone location, 
especially when located so close to the ground surface. At low frequencies the prediction shows an extra 
constructive and destructive pattern relative to the experiment. Overall, in the far-held at R/D = 69.23 and 
at heights high above the ground, the ground effects have been captured by the prediction method. 

The ground shadow distance from the nozzle exit to the microphone is set at R/D = 46.15 and the 
microphone heights are set at 3.69, 1.52, and 0.0061 m. Comparisons between predictions and experiment 
are shown in Fig. 6. At the heights of 3.69 m and 1.52 m in Fig. 6(a) and 6(b) it is shown that both the 
overall magnitude and peak frequency of the spectrum are within 2 to 4 dB of the measurement, except in 
the 1.52 m case where there is considerable discrepancy near the first instance of destructive interference. At 
higher frequencies the prediction is typically 0 to 6 dB lower than experiment. At the ground microphone 
position the predictions do not compare favorably. The peak frequencies of both predictions are higher than 
measurement by approximately a few hundred Hz and 2 to 6 dB lower. Once again a constructive interference 
pattern in the prediction is seen near 30 Hz. 

For the final predictions the microphone positions are moved well into the near-field of the jet. Near-field 
acoustic predictions from jet flows are often very difficult. The shadow distance from the nozzle exit to the 
microphones is set at R/D = 23.08. Figure 7(a) shows the predictions when the microphone height is at 
3.89 m. The peak magnitude of the impedance based model agrees very well with experiment and the hard- 
wall acoustic model is over-predicted by one dB. The peak frequency is approximately 450 Hz. The lowest 
frequency of the experiment is 40 Hz and is too high to observe any destructive interference. Unfortunately, 
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the high frequency fall-off of the predictions are larger than the measurement. 

A second near-field prediction is shown in Fig. 7(b) where predictions are compared with experiment at 
R/D = 23.08 m and a height of 1.52 m. The predictions are nearly identical except for the slight attenuation 
due to the ground. Both predictions over-predict the lower frequencies, are within 1 dB of the peak of 
the experiment, and greatly under-predict the energy above 900 Hz. In the experiment there is very little 
apparent constructive or destructive interference. A higher fidelity convective wave equation Green’s function 
solver could potentially correct the predictions in terms of the increased fall-off at higher frequencies and 
also correct for the extra interference at low frequencies. 

A final prediction is shown in Fig. 7(c). The microphone height is 0.0061 m and the shadow distance 
is R/D = 23.08 m. This is a difficult position to make an accurate prediction as the observer is on the 
ground and is within the near-field of the jet. The peak frequency of the hard-wall prediction is 100 Hz 
higher than the experiment and the impedance peak frequency is 200 Hz lower than the experiment. The 
attenuation between the hard-wall prediction and the impedance prediction is much larger than previously 
seen. Unfortunately, the prediction using a ground surface with realistic impedance greatly under-predicts the 
measurement. It should be noted how extremely sensitive to the microphone position and ground impedance 
the predictions are. Even a slight change of microphone height of 0.001 m near the ground can greatly 
influence the constructive-destructive interference patterns. Small changes in the impedance can also greatly 
change the predicted values. In the near-field of the jet, roughness or small irregularities on the ground 
surface can greatly change the experimental result. This causes obvious difficulty when predictions are 
assumed about a perfectly flat surface and are amplified in the near-field of the jet. Overall, the predictions 
are satisfactory considering the sensitivity of the problem, measurement, and conducted in the near-field at 
a jet emission angle that is relatively less understood than the sideline. 

Conclusion 

An acoustic analogy is formulated and the arguments are the vector Green’s function of the LEE and 
equivalent sources. The latter involves the two-point cross-correlation of the turbulent velocity fluctuations 
and scaling terms based on a steady RANS solution. The vector Green’s function of the LEE is written in 
relation to Lilley’s 4 equation, which is approximated with matched asymptotic solutions about the Green’s 
function of the convective Helmholtz equation. The Green’s function of the Helmholtz equation is found 
analytically for the infinite flat plane. Large savings have been made in terms of computational cost with a 
minimal loss of accuracy relative to finding the full solution of the vector Green’s function of the LEE. The 
new model has been exercised and compared with the experiment of Miles.' The predictions are shown to 
capture jet noise ground effects that are characterized by constructive and destructive interference patterns 
in the mid- and far-field and capture overall trends in the near-fielcl. However, there are discrepancies with 
respect to measurement of the magnitudes and frequencies of select constructive and destructive interference 
patterns. 

The choice to study ground effects is due to the availability of experimental data and the ease of finding 
the appropriate Green’s function. This model formulation allows a multitude of well-known Green’s function 
solvers to account for scattering of jet mixing noise as long as the plume is not highly influenced by the 
scattering surface. The problem has been greatly simplified by relating the vector Green’s function of the 
LEE to the Green’s function of the convective Helmholtz equation while still retaining the benefits of Lilley’s 
equation that accounts for jet shear layer refraction and amplification. 
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Figure 1. The coordinate system showing the nozzle, the ground plane, and the observer. The positive y-axis 
points away from the ground below the nozzle. The origin of the coordinate system resides at the center of 
the nozzle exit. 



x/D 


Figure 2. A slice of the computational domain showing the convergent contoured nozzle. The x- and y-axis 
have been normalized by the nozzle exit diameter. The nozzle geometry corresponds to the NASA Glenn 
Research Center SMCOOO nozzle. 
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x/D 

Figure 3. Contour maps of the steady RANS solution closed by the Menter 31 SST turbulence model. The 
coordinates are normalized by the nozzle diameter. In the top half plane contours of the streamwise velocity 
component, u m/s, are shown. In the bottom half plane contours of turbulent kinetic energy, TKE (m 2 /s 2 ), 
are shown. 



Figure 4. Free-held predictions of the SAE ARP876 35 method contained within the Aircraft Noise Prediction 
Program (ANOPP) and the developed acoustic analogy. The jet conditions are Mj = 0.95, TTR = 1.00, D = 0.3302 
m, R/D = 100, 'ijj = 145 degrees (inlet angle), and azimuthal angle <j> = 90 degrees. 
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(a) (b) 



(c) 


Figure 5. Ground shadow distance R/D = 69.2308 and microphone heights a) 3.89 m, b) 1.52 m, and c) 0.0061 
m (ground). 
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(a) (b) 



(c) 


Figure 6. Ground shadow distance R/D = 46.1538 and microphone heights a) 3.89 m, b) 1.52 m, and c) 0.0061 
m (ground). 
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(a) (b) 



(c) 


Figure 7. Ground shadow distance R/D = 23.0769 and microphone heights a) 3.89 m, b) 1.52 m, and c) 0.0061 
m (ground). 
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